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I.  INTRODUCTION 


If  one  measures  the  amplitude  of  the  wave  reflected  from  a  junction 
placed  in  double  hole  SQUID  fashion  across  a  waveguide,  with  provisions  to 
voltage  bias  the  junction  (Resistive  SQUID),  then  the  in  phase  and  out  of  phase 
components  vary  with  bias  as  shown  in  Fig.  1.  The  RSJ  model  qualitatively 
predicts  an  impedance  variation  of  this  type.  However,  in  previous  analyses 
the  out  of  phase  component  (or  reactive  portion)  appeared  only  when  the  junc¬ 
tion  oscillations  were  locked  to  the  probing  signal.  Reported  impedance  mea- 
2 

surements  as  well  as  our  own  observations  on  point  contacts  did  not  show 
any  evidence  of  locking.  As  a  matter  of  fact,  in  our  experiments  the  self  oscil¬ 
lations  could  occasionally  be  observed  as  hash  on  the  center  portion  of  the 
structure  in  Fig.  1.  Thus  locking  definitely  did  not  occur  and  another  inter¬ 
pretation  of  the  observed  impedance  variation  with  bias  is  required.  A  sim¬ 
ple  explanation  may  be  provided  by  the  RSJ  model  if  terminations  at  i  ’  er 
frequencies,  particularly  that  at  the  difference  frequency  between  applied 
signal  and  self -oscillations,  are  taken  into  account.  This  interpretation  has 
been  explored  by  a  small  signal  analysis,  which  is  described  in  this  com¬ 
munication.  It  reveals  that,  by  measurement  of  the  impedance  of  a  biased 
junction,  a  number  of  junction  and  circuit  parameters  could  be  determined, 
which  otherwise  would  not  be  accessible.  It  appears  that  these  measurements, 
together  with  an  I-V  characteristic,  would  serve  to  completely  characterize 
a  junction  and  its  immediate  embedding  network.  The  reason  for  this  is  the 
fact  that  with  the  junction  self -oscillations  we  are  able  to  provide  a  frequency 
sweep  from  essentially  zero  frequency  to  frequencies  very  large  compared 

* 

M.  Kaminaga,  "rf  Impedance  of  Superconducting  Weak  Links,  "  J.  Phys.  Soc. 
Japan,  to  be  published. 

*F.  Auraches  and  T.  Van  Duzer,  "rf  Impedance  of  Superconducting  Weak 
Links,  "  JAP  44,  848  (1973). 

^D.  E.  Claridge,  "rf  Current  Dependence  of  8.9  GHz  Impedance  of  Tantalum 
Point-Contact  Josephson  Junctions",  JAP  48,  1762  (1977). 


VOLTAGE  BIAS  OR  JUNCTION  FREQUENCY 


Fig.  1.  Change  is  reflection  when  junction  bias 
is  swept  through  probing  frequency  sig¬ 
nal  at 


with  the  operating  frequency;  and,  of  course,  frequency  is  an  essential 
parameter  in  circuit  evaluation.  No  other  test  system  exists  with  the  fre¬ 
quency  versatility  provided  by  a  Josephson  junction.  Thus,  it  appears  that 
the  impedance  evaluation  could  lead  to  quantitative  determination  of  series 
inductance,  junction  capacitance,  junction  ohmic  resistance  at  operating 
frequency,  low  frequency  (DC)  loop  termination,  and  fundamental  pump  ampli¬ 
tude.  An  I-V  characteristic  would  furthermore  provide  critical  current  and 
junction  impedance  at  DC;  and  finally,  the  frequency  profile  of  the  self¬ 
oscillations  at  operating  frequency  would  permit  determination  of  the  noise 
power  in  the  circuit.  This  provides  a  rather  complete  junction  description, 
lacking  essentially  only  the  spectrum  of  higher  harmonics  in  the  circuit. 
However,  since  the  effects  of  subharmonic  frequencies  of  the  bias  oscillations 
are  observable  in  experiment,  a  more  extended  analysis  might  help  to  give 
some  information  on  other  harmonics  in  the  circuits  as  well. 

However,  there  exist  various  difficulties  with  the  analysis  presented 
here,  and  I  have  not  been  able  to  come  to  a  satisfactory  solution  of  these 
difficulties.  They  are  pointed  out  in  the  respective  sections  of  this  note.  As 
usual,  refuge  is  sought  in  computer  simulations,  which  should  be  better 
suited  to  treat  the  strong  nonlinearities  in  the  device.  These  studies 
are  presently  in  progress  on  a  computer  with  the  aim  to  find  out  the  validity 
of  the  results  developed  in  the  analysis.  In  case  the  numerical  results  bear 
out  the  analysis  we  should  have  some  good  junction -circuit  diagnostics  at 
hand;  otherwise,  the  exercise  might  serve  to  at  least  give  some  indication  of 
the  qualitative  functioning  of  the  circuit. 
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II.  ANALYSIS 


We  briefly  derive  the  small  signal  approximation  of  the  impedance  for 

3 

the  RSJ  model  as  it  has  previously  been  done  by  Vistavkin,  et  al.  and  then 
include  the  effect  of  the  terminations.  We  start  with  the  circuit  equation 


+  I  sin  <l>  =  I  +  I  sin  uu  t 
c  r  o  s 


(1) 


with  the  variables  explained  in  Fig.  2.  We  assume  I  «  Iq.  Iq  causes  self¬ 
oscillations  at  uuo  with  phase  pQ(t).  The  phase  is  modified  by  the  small  signal 
pj(t)  «  <pQ(t).  The  linear  1st  order  approximation  in  reduced  units  is 


+  cos  <pQ( t)*  «/>1  =  i  (2) 

where  i  =  I/I  and  time  is  measured  in  units  of  cu  *  =  9  /I  RT.  Standard  inte- 
c  ^  c  o  o  J 

gration  and  the  fact  that  cos  =  -0^/0^  =  -  ^  fn  <p'Q  leads  to  the  small  signal 
voltage  (in  units  I  RT) 

C  J 


V  =  1  +  f  Jo  dT  (3) 

The  second  term  always  gives  a  real  (resistive)  contribution  which  can  become 
negative.  Note  that  it  describes  a  double  mixing  process  with  the  self-oscil¬ 
lations  by  two  90°  phase  shifting  processes:  (a)  integration  and  (b)  multiplica¬ 
tion  with  a  90°  phase  shifted  voltage  contribution  (when  only  the  fundamental 


A.  N.  Vistavkin,  V.  N.  Gubankov,  L.  S.  Kuzmin,  K.  K.  Likharev, 

V.  V.  Migulin  and  B.  K.  Semenov,  "S-c-S  Junctions  as  Non-Linear  Elements 
of  Microwave  Receiving  Devices",  Rev.  Physique  and  Appliquee,  9_,  79 
(1974). 
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frequency  terms  are  considered).  The  fact  that  v,  and  thus  the  impedance 
v/i,  is  determined  by  the  large  voltage  variation  due  to  self  oscillation  and 
its  derivative  suggests  that  the  above  equation  is  of  more  general  validity 
than  applicable  only  for  sin  p  dependences.  In  fact  it  can  be  shown  that  it 
holds  for  any  periodic  I  vs.  p  dependence  if  an  accommodation  between  super 
and  normal  current  is  enforced  by  the  external  circuit  (See  Appendix). 

Using  Eq.  (3),  one  obtains  an  impedance  matrix  (in  units  of  Rj)  for  the 
3  main  idler  components  (Fig.  3) 


where  it  is  assumed  that  small  currents  are  injected  at  all  three  frequencies 

of  interest  and  only  terms  have  been  retained  which  become  large  for  <d  ^  = 

uj  -  cd  -•  0,  or  1st  order  terms  if  no  singular  term  occurs.  The  matrix  holds 
so 

for  the  analytically  tractable  RSJ  model  with  I  =  I  sin  p,  where  the  terms 
—  .  _  i~Z  S  C 

(a  -  v),with  a  =  Iq/Ic  and  v  =  ^a^-  1,  represent  the  Fourier  component  at  u.'  , 
the  pump  frequency.  In  what  follows,  we  replace  (a  -  v)  with  a  more  general 
amplitude  coefficient  y  =  — 5 — . 

Using  standard  methods  after  inclusion  of  terminations  (where  the 
voltage-current  relation  is  enforced  by  the  external  circuit)  one  obtains 
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from  which  the  impedance  of  the  circuit  at  oj  is  obtained  (with  a  =  =  =  — — 

s  v  R  _ 

where  is  the  slope  of  I-V  characteristic):  ^ 


We  now  specify  the  termination  at  U)_(as  Z_{  =  ju^L^  (in  units  Rj),  where 
Lj-jq  is  the  DC  inductance  of  the  circuit  closing  the  junction  current.  Assum¬ 
ing  for  now  Z  ^  ~  ®,we  obtain 


ju,-iLPC  _y_ 
°  +  ja,-lLDC  "c 


(7) 


This  is  the  junction  impedance  in  the  vicinity  of  (D  u)  .  For  L  J  =  |u>  -  a) 

o  s  1  - 1  'so 

>:>  Rn^Lnr  the  Parenthesis  vanishes  and  Z  Rx.  At  w  .  =  0  the  reactive 
i  Uy  S  J  .  -1 

term  is  -j  y  LDC^RD’  The  "half  width"  holds  for  |af  =  RD/LDC>  From 
both  these  relations  which  can  be  determined  experimentally,  L  /R  and  y 

JJU 

can  be  determined. 

The  above  quantities  can  readily  be  read  off  a  Smith  chart.  We  assume 
Zg  to  be  in  series  with  a  lead  inductance,  L  .  Thus,  for  w  »  R  /L  and 
*  0  in  Eq.  7,  a  point  (Lg,  Rj)  is  determined.  When  sweeping  wq  through  (ug> 
the  {  }  term  assumes  finite  values  and  a  circle  is  traced  as  shown  in  Fig.  4. 
Both  circle  diameter  in  terms  of  ouT  „/Z  ,  and  fou^V  -  JZ?)/Z  =  R^/L,^  in 
terms  of  junction  bias  voltage,  are  obtained. 


While  the  result  presented  is  in  qualitative  agreement  with  experimental 
observations,  the  derivation  contains  a  serious  flaw.  In  case  a  real  portion  is 
considered  in  the  impedance  Z  ^  of  Eq.  (6),  no  analytic  expansion  is  obtained 
since  a  term  tu  remains  in  the  denominator,  and,  consequently,  the  term 
grows  beyond  bounds.  Of  course,  in  reality  the  junction  generator  at  m  ^  has 
an  internal  impedance  preventing  singular  behavior.  One  might  suspect  that 
as  long  as  the  external  real  load  at  ui  j  remains  small  compared  with  R^.,  the 
analysis  might  give  a  reasonable  description.  Presently,  numerical  simula¬ 
tion  is  hoped  to  shed  some  light  on  this  problem. 
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III.  EXPERIMENTAL 


In  our  10  GHz,  broad  band,  Zq  5fi  waveguide,  we  obtain  a  reflection 
response  of  just  the  type  discussed  above  (see  Fig.  5).  Here  we  observe 
reflected  vectors  for  constant  input  power  at  .  Thus,  consider  the  reflected 
vector  normalized  and  the  Smith  Chart  evaluation  applies  exactly.  The  ob¬ 
served  inductance  change  is  indeed  negative.  However,  it  is  also  apparent 
that  the  measured  circle  in  response  to  bias  variation  is  somewhat  distorted. 

In  what  follows  we  indicate  some  reasons  for  this. 

A.  We  first  consider  the  termination  at  to  Ignoring  the  series  inductance, 
Lg,  the  junction  near  oug  is  essentially  terminated  by  a  short.  Deviations  to 
either  higher  or  lower  frequencies  will  result  in  capacitive  or  inductive  load¬ 
ing  of  the  junction  at  u>  We  have  simulated  the  short  with  loading  by  a  tank 

circuit,  introducing  at  id  a  load  of  form  Z  ^  =  j  Q  —  where  A  =  cu  n  -  au  = 

s  -C.  (DS  -Z  s 

2u)  j.  This  idler  loading  is  evaluated  with  Eq.t,(6)  and  results  in  a  circle 
distortion  as  shown  in  Fig.  6.  We  see  distortion  towards  larger  vectors  for 
with  no  effect  at  u)q  =  m  in  qualitative  agreement  with  experimental 
indication  (Fig.  5). 

B.  The  clearly  observable  distortions  near  the  origin  of  the  circle  (curl -in 
at  Z  =»  Lg,  R j)  are  not  described  by  the  model.  Since  our  approximation  holds 
only  for  u)Q  **  u;g,  this  is  not  particularly  surprising. 

C.  Since  we  are  dealing  here  with  a  "variable"  inductance  at  u)  it  should 

s 

enable  one  to  determine  the  capacitance  of  the  junction.  If  the  assumption  is 

made  that,  in  a  circuit  as  sketched  in  Fig.  7a,  the  capacitance  does  not  affect 

the  general  junction  response  with  bias  as  determined  by  the  termination  at 

u)_jibut  enters  only  into  the  external  circuit,  then  it  can  be  shown  by  simple 

circuit  analysis  that  the  response  circle  is  tilted  on  the  Smith  chart  as  shown 

in  Fig.  8.  Considering  only  the  circle  origin*  and  the  point  for  w  =  w  ,  we 

OS 

*  ”  ’  ; 

Circle  origin  is  that  point  obtained  for  the  ideal  model  where  w  -•  o  or 

u  »  u  .  ° 

o  s 


-17- 


fRBCSOIMO  PiOB  m.AKy 


can  now  read  off  four  values  (two  for  the  origin,  two  for  w  =  cj  ).  These 

o  s 

values  determine  real  and  imaginary  parts  of  the  expressions  derived  from 
the  circuit  of  Fig.  7a: 

2 


The  result  is  four  equations,  which  are  sufficient  to  determine  the  unknowns 
L*S>  Rj>  C  and  X  of  Fig.  7a.  From  X  and  the  measured  half  width  we  deduce 
Rd/ldc  and  the  effective  pump  amplitude.  Clearly,  if  one  can  also  plot  the 
I-V  characteristic,  a  complete  characterization  of  the  junction  results. 

Thus  the  small  signal  analysis  based  on  the  RSJ  model  gives  a  descrip¬ 
tion  in  fair  qualitative  agreement  with  experiment.  The  difficulty  encountered 
with  the  assumption  in  C  above  is  touched  upon  in  the  appendix. 


(l  +  udCRj)2  +  (1  +  oiCX)2 


+  j  u>Lc 


(uCRj  +  X(i  +  ujCX) 

1  +  <dCRt)2  +  (1  +  ujCX)< 
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Fig.  8. 


Modification  of  idealized  response  (Fig.  4)  by 
capacitive  loading. 


APPENDIX 


Here  we  show  to  what  extent  a  periodic  supercurrent-phase  relation 
always  results  in  a  real  contribution  to  biased  junction  impedance.  Suppose 
the  junction  (with  no  capacitance)  is  driven  by  a  DC  current  in  RSJ  model 
fashion;  the  differential  equation  for  the  phase  evolution  then  is 

<t>  +  f(0  )  =  const  (IB) 

ro  ro 

when  written  in  reduced  units  and  when  the  usual  sin  0q  term  is  replaced  by 
an  arbitrary  periodic  function,  f(0Q).  The  response  to  a  small  ac  current 
injected  in  addition  to  the  DC  current  results  in 

<p'Q  +  0^  +  f(0Q  +  0^)  =  const  +  Ai  (2B) 

or,  after  expansion  to  first  order  and  use  of  Eq.  (IB), 

*'i  +  a?loh  =  Ai  (3B) 

On  the  other  hand  df/d0|Q  may  also  be  obtained  from  Eq.  (IB)  by  differentia¬ 
tion  with  respect  to  time: 


0"  +  0'  =0 
ro  d<f>  To 


(4B) 


from  which 


(5B) 
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and  consequently  Eq.  (3B)  may  be  written 


-(&  =  Ai 


Standard  methods  for  solving  a  first  order  linear  differential  equation 
<t>'l  +  g(t)p^  =  Ai  result  in 


'/g  dt  [  /g  dt 


<t>^)  =  -  exp  J  exp  A^  dt  +  Q  exp 


-fg  dt 


or,  with  Eq.  (6B ) 


dt  +  Q  <p' 


-  »i<‘>  -  <  ffo  dt  +  Ai 


where  we  choose  initial  conditions  such  that  Q  =  0.  Assuming  only  the  funda¬ 
mental  Fourier  components  of  9/  and  <f> "  to  be  of  significance,  it  is  easily 
verified  that  the  first  term  of  Eq.  (7B)  generates  in  phase  components  only. 
We  now  briefly  consider  the  more  general  case  where  the  junction  is  em¬ 
bedded  in  the  circuit  of  Fig.  7b.  The  starting  equations  are 


CV  +  17  +  =  h 


L  -37  I  +  I  R  +  V  +  V  =  0 
dt  T  To  o 
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After  integration  of  Eq.  (9B)  and  insertion  into  Eq.  (8B)  we  find 


cv  +  rt  +  + 


i  R  r  r 

L  exp  -  —  1  J  v  exP 


t  dt  = 


R 


(10B) 


We  approximate  the  4th  term  by  setting  in  it  V  =  v  +  v  cos  <u  t  and  obtain 


(11B) 


With  u)  L  «  R  we  may  neglect  the  cos  u)  t  contribution.  Then  the  second 

CJ  O  O 

portion  of  Eq.  (11B)  becomes 


uu  L 

o 


R 

o 


v  cos  a)  t  dt 
o 


Eq.  (10B)  can  then  be  written  in  reduced  units 


0)  L 

0C0,/  +  <f>'  +  S(0)  +  — 

o 


V  -  v 
o 

R  I 


o  c 


(12B) 


whe  re 

9  =  hC  /I  2e 

c  c 

0  =  h/ZeIcL 
Note  that 


a)  L 
o 

R 

o 


i 

9 


0)  R  _ 

o  J 

Cl)  R 

c  o 


(13A) 


-27- 


which  for  a  properly  loaded  junction  with  u0Lj  RQ  is  of  order  one;  therefore 

even  for  small  <<  Rq  the  term  is  usually  not  negligible.  Furthermore 

v  in  the  right  hand  side  changes  to  second  order  only  with  application  of  a 

small  disturbance  in  V  -•V  +  Av  cos  u)t.  If  we  set  8  =0,  the  linearized 

o  o  s  c 

small  signal  approximation  of  Eq.  (12B)  is  of  the  same  form  as  Eq.  (IB)  and 
we  may  again  write  the  solution 


v  =  0't(t) 


/,  Av 
i  _ s_ 

I  R 

o  c  o 


Av 


cos  (i)  t  dt  +  ■=■ 
s  I 


R 

c  o 


COS  0)  t 

s 


(13B) 


Clearly  the  same  conclusions  apply  as  before  for  the  current  driven  junction. 
The  junction  impedance  is  essentially  real  and  modification  due  to  idler  load¬ 
ing  as  described  in  the  main  test  should  apply.  The  L  in  series  with  the  junc¬ 
tion  does  not  modify  the  basic  result. 

Here,  the  difficulty  arises  in  asserting  that  the  analysis  applies  also 
when  O  or  6c  41  O.  When  only  supercurrent  and  an  ohmic  path  is  pro¬ 
vided  by  the  junction,  with  the  requirement  that  the  sum  total  remains  constant, 
the  above  analysis  applies.  Clearly,  the  reasoning  cannot  hold  any  longer  if 
an  additional,  reactive  bypass  is  provided.  Again,  since  the  mathematics 
become  complicated,  it  might  be  more  economical  to  test  with  numerical 
simulation.  The  results  of  this  simulation  will  be  reported  at  a  later  time. 
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LABORATORY  OPERATIONS 

The  Laboratory  Operations  of  The  Aerospace  Corporation  is  conducting 
experimental  and  theoretical  investigations  necessary  for  the  evaluation  and 
application  of  scientific  advances  to  new  military  concepts  and  systems.  Ver¬ 
satility  and  flexibility  have  been  developed  to  a  high  degree  by  the  laboratory 
personnel  in  dealing  with  the  many  problems  encountered  in  the  nation's  rapidly 
developing  space  and  missile  systems.  Expertise  in  the  latest  scientific  devel¬ 
opments  is  vital  to  the  accomplishment  of  tasks  related  to  these  problems.  The 
laboratories  that  contribute  to  this  research  are: 

Aerophysics  Laboratory:  Launch  and  reentry  aerodynamics,  heat  trans- 
fer<  reentry  physics,  chemical  kinetics,  structural  mechanics,  flight  dynamics, 
atmospheric  pollution,  and  high-power  gas  lasers. 

Chemistry  and  Physics  Laboratory:  Atmospheric  reactions  and  atmos- 
pheric  optics,  chemical  reactions  in  polluted  atmospheres,  chemical  reactions 
of  excited  species  in  rocket  plumes,  chemical  thermodynamics,  plasma  and 
laser-induced  reactions,  laser  chemistry,  propulsion  chemistry,  space  vacuum 
and  radiation  effects  on  materials,  lubrication  and  surface  phenomena,  photo¬ 
sensitive  materials  and  sensors,  high  precision  laser  ranging,  and  the  appli¬ 
cation  of  physics  and  chemistry  to  problems  of  law  enforcement  and  biomedicine. 

Electronics  Research  Laboratory:  Electromagnetic  theory,  devices,  and 
propagation  phenomena,  including  plasma  electromagnetics:  quantum  electronics, 
lasers,  and  electro-optics;  communication  sciences,  applied  electronics,  semi¬ 
conducting,  superconducting,  and  crystal  device  physics,  optical  and  acoustical 
imaging;  atmospheric  pollution;  millimeter  wave  and  far-infrared  technology. 

Materials  Sciences  Laboratory:  Development  of  new  materials;  metal 
matrix  composites  and  new  forms  of  carbon;  test  and  evaluation  of  graphite 
and  ceramics  in  reentry;  spacecraft  materials  and  electronic  components  in 
nuclear  weapons  environment;  application  of  fracture  mechanics  to  stress  cor¬ 
rosion  and  fatigue-induced  fractures  in  structural  metals. 

Space  Sciences  Laboratory:  Atmospheric  and  ionospheric  physics,  radia- 
tion  from  the  atmosphere,  density  and  composition  of  the  atmosphere,  aurorae 
and  airglow;  magnetospheric  physics,  cosmic  rays,  generation  and  propagation 
of  plasma  waves  in  the  magnetosphere;  solar  physics,  studies  of  solar  magnetic 
fields;  space  astronomy,  x-ray  astr  inomy;  the  effects  of  nuclear  explosions, 
magnetic  storms,  and  solar  activity  on  the  earth's  atmosphere,  ionosphere,  and 
magnetosphere;  the  effects  of  optical,  electromagnetic,  and  particulate  radia¬ 
tions  in  space  on  space  systems. 
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